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Atmospheric free fall with quadratic resistance when considered from the point of view of 
the strictly isospectral construction of supersymmetric quantum mechanics is a powerful 
analogy for the class of Darboux cosmological fluids that I introduced in a previous work. 
Here, using this analogy, it is shown that the possible accelerating cosmological effects 
are due to the quadratic negative friction properties with which the Darboux fluids are 
endowed through the mathematical construction. 

1 - It has recently been realized by Faraoni,0 that the equations describing 
barotropic FRW cosmologies can be combined in a simple Riccati equation leading 
in a much easier way to the textbook solutions for the FRW scale factors. In a 
previous work, using mathematical schemes belonging to supersymmetric quantum 
mechanics, I have employed the FRW Riccati equation to introduce a so-called 
Darboux class of cosmological fluids with time-dependent adiabatic indices 7(77; A)J3 
This type of fluids could provide an explanation for the recent data from high- 
redshift Type la (hydrogen-deficient) supemovae indicating that the universe has 
now an accelerating large scale expansion.13 In this note, I would like to clarify 
the nature of the effects produced by the mathematical scheme through the close 
analogy with the atmospheric resistive free fall. 

2 - Faraoni obtained the following cosmological Riccati equation 

— +cu 2 + kc = Q (1) 

for the log derivative of the FRW scale factor u = da ^ dri . The independent variable 
is the conformal time 77, c = §7 — 1 is assumed positive, and n — 0, ±1 is the 
curvature index of the flat, closed, open FRW universe, respectively. 

From the mathematical point of view, the FRW Riccati equation is of constant 
coefficients and therefore is directly integrable. In general, for a Riccati equation of 
the type 

^-=ay 2 +0, (2) 
ax 
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the solutions can be written as follows 



(i) for | > 



l/ = y£tan (J^px + ^y or y = y ^ cot ^- a/^^ + 4>i j t^' 



(ii) for | < 

-^tanh (J~Px + (p^j , or y = Jj^coth^J~Px + tp^ ■ ( 4 ) 
One can see that the cosmological cases correspond to 

a 1 

a = -c, /3 = -kc ==>-- = - (5) 

and therefore the two mathematical cases are directly connected to the type of FRW 
universe. The case k = is presented in H|. 

The time-dependent barotropic indices 7^(77; A) have been introduced through 
the so-called strictly isospectral Darboux technique (SIDT). The latter can be inter- 
preted as a three-step procedure in which one goes from a physical Riccati problem 
to an unphysical one by changing the free term and the sign of the derivative term 
and then comes back to another physical problem similar to the original one by 
the same type of changes. However, in the third-step problem, unphysical effects of 
the intermediate problem are maintained as a consequence of the applied scheme, 
essentially because the general Riccati solution of the intermediate problem is used 
in the last step as well. The very illustrative example of the classical Newton motion 
with quadratic resistance is presented below since it both helps understanding the 
three-step character of SIDT and provides further insight in the nature of Darboux 
cosmological fluids. 



3 - Newton's law of force for classical motions with quadratic resistance, such as 
gravitational free fall in terrestrial atmosphere and similar motion of electric charges 
in constant electric fields, is Riccati's equation without linear term. About 40 years 
ago, the free fall in quadratic resistive atmosphere has been studied in some detail 
by Davies in his book.El The equation of motion is of the following type 



dt 2 



mg 



(6) 



which can be written 



dv 



ev 



(7) 
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where e = — is the friction coefficient per mass unit. Again, comparing this New- 
tonian case with Eq. (2) one gets 

a — 6 

a = -e, = g =► - = — < 0, (8) 

P 9 

and therefore the solution is given by 

v(t) = ^tanh {y/qjt + if!) , (9) 

which is Davis' result. The falling velocity v(t) presents a substantial time de- 
pendence - acceleration stage - only for a short period and rapidly turns into the 
limiting constant velocity V = yp£. The phase ipi is related to the initial velocity. 
It is not important for the general discussion here and therefore we take <pi = 0. 

Let us sketch now shortly, the three steps of SIDT for the resistive atmospheric 
free fall: 

(i) The initial equation is ^ = g — ev 2 representing the normal resistive falling 
motion. 

(ii) One shifts to the supersymmetric 'unphysical' partner motion of the form 

^ = 92(f) - ev 2 , (10) 

where 

g 2 (t) =.g[-l + 2tanh 2 (V^t)] . (11) 

The negative sign in Eq. (10) shows that the motion is upwards and through a 
fluid medium endowed with friction of negative friction coefficient (i.e., antifriction) 
as the driving force, plus another force g% (t) of driving character for t < T = 
— |=Arctanh-^ and turning friction-like afterwards. 

The first two steps are connected to each other through the particular solution 
of the initial normal motion as given by Eq. (9). 

(iii) The third step is the return back to an equation similar to the initial one, 
namely 

(11 7 

-=g 3 (t;X)-ev 2 , (12) 
where the time-dependent forces gs(t; A) are given by 



93 (t; A) =g 



(13) 



where Ib(t) — cosh 2 xdx. 

The second and the third steps are also connected to each other through the 
general solution of the 'unphysical' problem. Thus, although the motion given by 
Eq. (12) looks normal, its solution v(t) corresponds to a rising on the vertical rather 
than a fall. Exactly as in the falling case, the rising has a short accelerating stage 



4 Haret C. Rosu 



followed by a uniform rising with constant limiting velocity. Indeed, for sufficiently 
large times, the driving antifriction force is at its maximum value —eV 2 , whereas 
giit — ► oo) = +g is a constant friction force and one gets the value of V from the 
balance of the two forces. 



4 - Eqs. (5) and (8) give the translation rules for the analogy between the cos- 
mological case and the atmospheric friction case. From Eq. (1) one sees that the 
normal cosmological friction force is cu 2 . The equivalent of the upward velocity is 
Ufi = a (rj-X) da "j^' > ^ , whereas the ones for g 3 (t; A) are the time-dependent indices 
7 K (77; A). The analogy also predicts that the accelerating period is limited in time. 
In this 'quadratic friction' picture, the cosmological SIDT as worked out in our 
main paper ^ is of unphysical nature as far as the 'antifriction' properties of the 
Darboux fluids should be invoked through the second step of SIDT. However, cos- 
mology is perhaps the only branch of science where unphysical effects are manifest. 
The aforementioned SIDT atmospheric rising looks completely unphysical but the 
expansion of the universe is a reality. 
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ERRATA AND ADDENDA 



DARBOUX CLASS OF COSMOLOGICAL FLUIDS WITH 
TIME-DEPENDENT ADIABATIC INDICES 



H.C. Rosu 

[Mod. Phys. Lett. A, Vol. 15, No. 15 (2000) 979-990] 



A number of corrections to the work in which I introduced one-parameter Darboux su- 
persymmetric barotropic fluids in cosmology are performed. As addenda, the flat FRW 
case, which was left apart, is presented in detail in the supersymmetric approach. 



Errata. 

On p. 981 

(i) In the second line of the last paragraph the tanh function should be coth. 
On p. 982 

(ii) The function w in Equation (7d) should carry the subindex f to emphasize 
that the solution Wf is not the same as in Eq. (7c), i.e., 

(d~77 ~ CUp ) (dr/ + CUp ) Wf = Wf _ c (~ u p + cu p) Wf = w £ _ c ' c K,susyWf = , 

{Id) 

(iii) In the unnumbered equation following (7d) the tan function should be 
squared, whereas the tanh should be coth squared. 

(iv) At the end of the same paragraph, the following sentence should be included: 
The solutions Wf are Wf = — , c , ,., and Wf = £ — 7 for k = 1 and k = — 1, 

1 1 cos(c7/-fd) 1 smn(c7?) 

respectively. 

(v) The general Riccati solution should be written %(ry) = U v {rj) — This 
leads to a free term equal to unity in the Bernoulli equation (10) and to the dropping 
of c everywhere in the argument of the integral T. In addition, in Eq. (11) c in the 
argument of should be dropped. 

(vi) In Eq. (12b), T CKi {j]) should be replaced by T K {j]). 

(vii) The left hand side of Eq. (13a) should be —kc k (t); A). 

On p. 990, the page number of Ref. 8 is 4850 and not 2850, whereas the e-print at 



Ref. 15 should be |hep-ph/0002297 
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Addenda. 

We would like to present shortly the flat case k = in the supersymmetric approach. 
For this case the analog of Eq. (4) at page 980 reads 



, 2 

a /a 



a +C UJ = °' (1) 
where c = 57 — 1. Writing b = |, Eq. (1) is rewritten in the form 

b + (c + l)b 2 = . (2) 

There are two cases. 

(i) For c + 1 7^ 0, the solution is b p = j4rr^, or in terms of the scale factor 

1 2 

a = aot c +! = aot 3 T , which is a standard result. 

(ii) For c+1 = (7 = 0), one gets b = 0, i.e., b p = | = const =H (in cosmology). 
Thus, we get the standard result a = ao exp(Ht). This case is completely degenerate 
from the point of view of the supersymmetric constructions. 

For the first case, the solutions b p enter factorizations of the following type 

(dT + (c + 1)bp ) (ir ~ (c + 1)bp ) w = * _ (c + 1)[bp + (c + 1)b p ]w = • (3) 

Thus, because of Eq. (2), we get w = 0. The solution is a linear function Wb oc t. 
The supersymmetric partner (or fermionic) equation with respect to Eq. (3) will be 

(di ~ (c + 1)bp ) (di + (c + 1)bp ) Wf = % ~ (c + 1)[ ~ bp + (c + 1)bp]wf = ■ 

., (4) 

This leads to Wf — ^Wf = 0. The linearly independent solutions are Wfi = t 2 and 

Wf 2 = t _1 . 

The general Riccati solution comes out from 

(c+1)b ^^ = (c+1)b; _^, 

whose general solution can be written down as b g (t) = b p (t) — ^j^j, where v(t) 
is an unknown function. Using this ansatz, one obtains for the function v(t) the 
following Bernoulli equation 

^ + 2(c+l)v(t)b p (t) = l, (5) 

that has the solution 

T (t) + A 
wg(t) 

where X (t) = J* w^(y) dy oc t 3 /3, and we consider A as a positive integration con- 
stant that is employed as a free parameter. 
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Thus, b g (t) can be written as follows 



b g (t;A)=b p (t) 



1 d 



c + ldt 



ln(J (t) + A) 



(7) 



d 
dt 



In 



Wb(t) 

[X (t) + A] 



(8) 



One can get a one-parameter family of quantities Pb(t; A) as follows 



P b (t;A) = (c+l)b 2 (t;A) + 



lbg(t;A) 



dt 



(9) 



2 d 2 
c + 1 dt 2 



ln(J (t) + A) 



(10) 



4w b (t)w b (t) 



+ 



2w£(t) 



(c + l)(J (t) + A) (c+l)(Jo(t) + A) 2 



(11) 



The interpretation of Pb is similar to that in the original paper, namely as the 
general Riccati contribution to any constant barotropic index 7. The corresponding 
parametric scale factors of the universe will be 



a b (t, A) 



Wb(t) 



[2b(t) + A] 



(12) 



Plots of 7(t; A) and a b (t; A) for the dust universe and for a particular value of A 
are presented in Figs. 1 and 2, respectively, together with the standard counterparts. 
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gamma(t; 0.9) 




Fig. 1 Time-dependent adiabatic index of a Darboux barotropic fluid of dust type for 
A = 0.9. There is a —1/t 2 singularity at the origin, after which the index becomes rapidly 

7 = 1. 




Fig. 2 Scale factors of a 7 = 1 barotropic fiat universe and of the corresponding 
one-parameter supersymmetric flat universe (for A = 0.9). 
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A one-parameter family of time dependent adiabatic indices is introduced for any given 
type of cosmological fluid of constant adiabatic index by a mathematical method belong- 
ing to the class of Darboux transformations. The procedure works for zero cosmological 
constant at the price of introducing a new constant parameter related to the time de- 
pendence of the adiabatic index. These fluids can be the real cosmological fluids that 
are encountered at cosmological scales and they could be used as a simple and efficient 
explanation for the recent experimental findings regarding the present day accelerating 
universe. In addition, new types of cosmological scale factors, corresponding to these 
fluids, are presented. 

1 - Recent data from type la supernovae that can still be considered as meager 
provoked tremendous interest since they indicate that the universe has now an 
accelerating large scale expansion,El (p + 3p < in Eq. (1) below). This supernovae 
revolution stimulated model builders to supply a remarkable amount of paper work 
on this new topic, in which such impressive words as quintessence can be found. 
On the other hand, simple explanations should always be taken into account on the 
background of ever more complicated modeling. In the following, I would like to 
act according to the latter rule and present a rather nontrivial contribution that 
can seemingly be related to accelerating universes. 



2 - The scale factor a(t) of a FLRW metric is a function of the comoving time 
t obeying the Einstein- Friedmann dynamical equations (EFDEs): 

a, 47rG. „ , 

- = ^-(P + 3 P ) (1) 



* This essay received an "honorable mention" in the Annual Essay Competition of the Gravity 
Research Foundation for the year 2000 - Ed. 
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k\ &TtGp K 



a 

where p and p are the energy density and the pressure, respectively, of the perfect 
fluid of which a classical universe is usually assumed to be made of, and k = 0, ±1 is 
the curvature index of the flat, closed, open universe, respectively. To solve EFDEs 
for a(t) one needs an equation of state for the fluid, i.e., a relationship between p 
and p. Very used ones are the barotropic equations of state (BES) 

P = (7 - 1)A (3) 

where 7 is the constant adiabatic index. As is well known, all the important cos- 
mological fluids assumed to be dominant along various cosmological epochs belong 
to the barotropic class. It is straightforward to get a(t) by direct integration for 
the flat case,_however there are some difficulties in the other two cases. Landau 
and Lifschitz a gave a general procedure to perform the integration, which is valid 
for any 7. On the other hand, forthe purposes of the present work, an equivalent 
procedure introduced by Faraoni U leading practically to the same results (however 
with key constraints) is very useful. 

3 - The alternative method of Faraoni to derive the scale factor of the universe 
is based on the Riccati equation and goes as follows. First, one combines EFDEs 
and BES yielding 

a, /a\ 2 ck 

-+c - +-=0, 4 
a \ a / & z 

where c = 57 — 1. The case k = is directly integrable i and I shall not pay 
attention to it in the present work though it can be included in the scheme. Eq. (4) 
is rewritten in the conformal time variable rj in the form 

— + (c-l)f— ) +ck = 0. (5) 
a V a I 

One can immediately see that by means of the change of function u = — the 
following Riccati equation is obtained (c > henceforth) 

u' + cu 2 + kc = . (6) 

Furthermore, employing u = - — one gets the very simple second order differ- 
ential equation 

w" + kc 2 w = . (7a) 

For k = 1 the solution of the latter is wi = Wi cos(c?7 + d) , where d is an arbitrary 
phase, implying 

aifo) = Ailcos^ + d)] 1 / , (8) 
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whereas for k = — 1 one gets w_i = W_isinh(c?y) and therefore 

a_i(r?) = A_ x [shivery)] 1 / , (9) 

where W±i and A±i are amplitude parameters. Eqs. (8) and (9) are the same so- 
lutions as in the standard procedure. According to Faraoni, the alternative method 
requires 7 be a constant. However, as I will show in the following there is still some 
room in this method for time-dependent adiabatic indices. 



4-1 will now apply the so-called strictly isospectral Darboux technique (SIDT) 
to Eq. (7a). This is a mathematical approach related to the factorizations of the 
second order linear differential equations that underlined supersymmetric quan- 
tum mechanics and the dynamical symmetry breaking ideas of Witten almost two 
decades ago.u There is a slight but important difference between Witten's approach 
and SIDT. While Witten's factorizations rely on the particular Riccati solution, the 
SIDT uses the general Riccati solution. This was first discussed by Mielnik in 1984. 
B The difference with respectto Mielnik is that I will employ the SIDT scheme at 
fixed zero 'energy' (see also 0). L refer the reader to my short survey of Darboux 
transformations for more detailsfl During the last few years, I have acquired some 
experience with SIDT extending its range of applications beyond supersymmet- 



ric quantum mechanics, to subject s s uch as Newtonian damping, U Fokker-Planck 
equation,!] and quantum cosmology.E3 It is my aim here to apply this mathematical 
scheme to the dynamics of the scale factor of a classical universe. 

First of all, despite c looks like a frequency for the function w, I shall interpret ±c 
as Schrodinger 'potentials'. This is preferable to ±c 2 for algebraic reasons coming 
out from the bunch of formulas that I present next. (Notice that we have a constant 
potential having a coupling constant equal to the potential itself!). Thus, Eq. (7a) 
will be considered as a Schrodinger equation at zero 'energy' (or more rigorously at 
zero factorization constant). 

The point now is that the Riccati solution u p = - — mentioned above is only the 
particular solution, i.e., u Pl i = — tan cry and u Pi _i = cothcry for k = ±1, respectively. 
The particular Riccati solutions are closely related to the common factorizations of 
the Schrodinger equation. Indeed, Eq. (7a) can be written 

w — c(— kc)w = (7b) 

and using Eq. (6) one gets 

d \ / d 



dry ' CUp ) ( dry ~ C% ) W = W ~ ° (Up + ° Up)W = ° ' (?C) 
To fix the ideas, we shall call Eq. (7c) the bosonic equation. On the other hand, 



tin my previous works, I have also used the acronym DDGR for double Darboux general Riccati 
to denote the same scheme.) 
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the supersymmetric partner (or fermionic) equation of Eq. (7c) will be 



— cu t 



Thus, one can write 



-+cu p ,w f 



W f - C(-U p + CUp)wf — W f - C • C re ,susyWf = . 



(7d) 



--K.susy 



(V) 



-u p + cu p 



c(l + 2tan 2 c?7) if n > 1 
c(-l + 2coth 2 cr?) ifft<l 



for the supersymmetric partner adiabatic index. The solutions Wf are Wf 



and Wf = 



for k = 1 and k = — 1, respectively. 



cos(c?7+d) 



sinh(c?7) 

To find the general Riccati solution, one is led to consider the following Riccati 
equation cu 2 — = cu 2 — whose general solution can be written down as 

Ugiv) — u p( 7 ?) — c vt^) ' wnere v (v) i s an unknown function. Using this ansatz, one 
obtains for the function v(?7) the following Bernoulli equation 



dry 



+ 2cv(t7)u p (?7) = 1, 



that has the solution 



v(r?) 



(10) 



(11) 



where X K (77) = f£ w 2 (y)dy, if we think of a half line problem for which A is a 
positive integration constant thereby considered as a free SIDT parameter. 
Thus, u g (7y) can be written as follows 



u g (f?;A) = Up (77) 



_d_ 

dr] 



In 



1 d 
c CI77 



ln(X K (77) + A) 



(12a) 



(126) 



The important result provided by the SIDT is the one-parameter family of adiabatic 
indices c K (r?; A) 



-KC K (ry; A) = cu 2 (?7; A) + 



dugfa A) 
d?7 



2 d 2 
c ar] z 



(13a) 
(136) 

c(X K (r,) + X) ' c(l K ( V ) + Xf ■ (13C) 
One can easily infer the following formula for the parametric adiabatic indices 

7«(??;A) 



= — KC — 



ln(Z K (77) + A) 
4w K (?7)w^(77) 2w4(7y) 



7«(»?;A) =7 + 



3k(3 7 - 2) d?] 2 



ln(J K (77) + A) 



(13d) 
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c«(?7;A) = c K , susy (r/) + --r^ln ( —. ^ ), ( 15 ) 



which I used for plotting. All c K (?7; A) have the same supersymmetric partner index 
c K , S usy(^) obtained by deleting the zero mode solution w K . They may be considered 
as intermediates between the initial constant index kc and the supersymmetric 
partner index c KjSUSy (?7). From Eq. (126) one can infer the new parametric 'zero 
mode' solutions of the universe for the family of barotropic indices c K (?y; A) as follows 

A so-called double Darboux feature of the SIDT can be inferred by writing the 
parametric family in terms of their unique "fermionic" partner index 

2d 2 / 1 
cdrf \w K (r];X) 

which shows that the SIDT is of the inverse Darboux type0 This is very impor- 
tant because it allows a two-step (double Darboux) interpretation, namely, in the 
first step one goes to the fermionic system and in the second step one returns to 
a deformed bosonic system. In addition, this interpretation helps me to introduce 
another conjecture, on which I commented in a different context elsewhere.t3 In par- 
ticle physics supersymmetric ideas forces one to introduce supersymmetric partners 
for each known particle. However, if one promotes SIDT at the level of a matter of 
principles, one can think of a scenario in which the supersymmetric partner system 
is unstable/unphysical and only the systems generated by means of the general Ric- 
cati solution are encountered in nature. Thus, in the cosmological context, I claim 
that only the time-dependent barotropic fluids of indices c K (r?; A) make the real ma- 
terial content of the universe. At the cosmological scale, the fluids of constant 7 
occur in the limit A — > 00 and therefore they look more as academic cases in this 
context. SIDT introduces a new time scale in zero-A cosmology, which is given by 
the integration constant A. In ordinary supersymmetric quantum mechanics, this 
constant is related to the contribution of the irregular zero mode entering the SIDT 
parametric zero modes. Extrapolating this fact to cosmology, one can say that this 
time constant gives the contribution of the 'unphysical' solution to the parametric 
modes given by eq. (14). In the common approach, the 'unphysical' solution is 
discarded by imposing initial conditions. 

Another comment I would like to make is that from the SIDT standpoint I have 
worked at both fixed zero factorization 'energy' and fixed coupling constant c of 
the Schrodinger 'potential'. Thus, this is different from both standard quantum 
mechanical framework, working at fixed coupling constant but variable energy, and 
the Sturmian approach, which works at fixed energy but variable coupling constant 
(depending only on quantum numbers) . From the strictly technical point of view, 
what I have introduced here is a supersymmetric class of solutions for the scale 
factor of the universe corresponding to a family of cosmological fluids connected to 
any given fluid of constant adiabatic index by the SIDT supersymmetric procedure 
at zero factorization constant. One can interpret the denominator in Eq. (14) as an 
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amplitude damping for these supcrsymmctric scale factors whose origin is the time 
dependence of the adiabatic index. Some plots for the matter-dominated universe 
(c = i) and radiation- filled universe (c = 1) are presented in Figs. l(a)-(d) and 
2(a)-(d), respectively. 

5 - In conclusion, a supersymmetric class of cosmological fluids has been intro- 
duced in this work possessing time-dependent adiabatic indices. Such fluids may 
provide a simple explanation for a currently accelerating universe and therefore 
they should be of considerable interest. In future modeling along the lines of this 
work the parameter A can be either a constant fixed e.g., by some variational prin- 
ciples or can be made dependent on other astroparticle parameters. At the general 
level, a richer dynamics of the scale factor of the universe is introduced by simple 
mathematical means without emphasizing the phcnomcnological features, which are 
always of considerable speculative origin. Moreover, the methods of this work can 
be easily applied to the majority of inflationary scenarios. 

I would also like to point that in supersymmetry breaking terms SIDT is p£ 
the unbroken type. Since several years, Witten questions this fundamental issueJ13 
Most recently, he remarked £il "In the standard framework of low energy physics, 
there appears to be no natural explanation for vanishing or extreme smallness of the 
vacuum energy, while on the other hand it is very difficult to modify this framework 
in a sensible way. In seeking to resolve this problem, one naturally wonders if the 
real world can somehow be interpreted in terms of a vacuum state with unbroken 
supersymmetry" . In the cosmological context of this work, there appears to be a 
simple technical explanation for the vanishing of the vacuum energy, because only 
in this case Eq. (4) reduces to a Riccati equation wherefrom one can proceed with 
SIDT, which is naturally of unbroken type and moreover it might be used to explain 
the supernovae la data as I argued herein. In the vacuum case, the supersymmetric 
one-parameter adiabatic indices are 7 K ,vac(?7; A) = — j^£^z ln(2' K (c7;) + A) , leading 
to plots similar to the other parametric barotropic fluids (See Figs. 3(a) and (b)). 

Finally, I mention that in the framework of supersymmetric quantum mechanics 
I published a direct iteration of this type of SIDT leading to multiple parameter 
'zero modes' cB, which can be adapted to cosmology if the approach presented here 
is proven relevant. 
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Fig. la. Zero mode scale factor of matter dominated universe for n — 1 and A = 0.9 
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Fig. lb. Zero mode scale lactor ol matter dominated universe for k — — 1 and 

A = 0.9 . 




Fig. lc. 7k(t?; A) for a matter dominated universe of k = +1 and A = 0.9 
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Fig. Id. 7k (??; A) for a matter dominated universe of n = 



-1 and A = 0.9 
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Fig. 2a. Zero mode scale factor of a radiation-filled universe for k = +1 and A = 0.6 
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2b. Zero mode scale factor of a radiation-filled universe for k = 



-1 and A = 
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Fig. 2c. 7k(?7; A) for a radiation dominated universe of k = +1 and A = 0.6 
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Fig. 2d. 7 re (rj; A) for a radiation dominated universe of n = — 1 and A = 0.6 . 
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Fig. 3a. 7k (?7; A) for a vacuum-dominated universe of n — 1 and A = 0.9 



Darboux cosmological fluids 




